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This paper describes an equalization procedure for systems using 
manually adjustable bump equalizers that is based on a mean-squared 
error criterion. We show that, in accordance with a Gauss-Seidel iteration 
process, gain adjustment always converges to the optimal value at which 
the minimum MSE of the equalized channel is obtained. Both zero 
forcing and MSE algorithms based on the Gauss-Seidel iteration method 
are derived, and hardware implementation of these algorithms is dis- 
cussed. According to the error reduction analysis, an equalizer composed 
of orthogonal networks requires only one iteration to bring the equalizer 
to the optimum state. For the bump equalizers used in the latest L5 Coaxial 
Carrier Transmission System whose Bode networks are semi-orthogonal, 
two to three iterations are shown to be sufficient to achieve the optimum 
gain settings in the mean-squared error sense. 

I. INTRODUCTION 

In this paper, a manual adjustment of bump equalizers is described 
which uses a mean-squared error (MSE) criterion. In the existing L4 
Coaxial Transmission System, the bump equalizers (realized with Bode 
equalizer networks 1 ) are used for line equalization and adjusted ac- 
cording to a zero forcing (ZF) algorithm. 2 This method results in an 
optimum equalization in the MSE sense, but only under certain very 
restrictive conditions with respect to the transmission response of the 
channel. The latest L5 Coaxial Transmission System, which provides 
up to 10,800 toll-grade long-haul message channels on a pair of 0.375- 
inch coaxial cables over 4000 miles, also includes bump equalizers for 
equalization of the 65-MHz bandwidth channel. The equalization 
method used in the L5 Coaxial Carrier Transmission System, however, 
minimizes the MSE of the channel deviation. It is found in practice 
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that the MSE algorithm has given a better equalization result than 
the ZF algorithm. 

Since a coaxial cable system shows relatively stable channel charac- 
teristic, the bulk of the L5 line equalization is accomplished by 
manually adjustable equalizers. Normally, the time-varying channel 
deviation in a cable system is mostly the result of seasonal tempera- 
ture variation and aging of the components in the system. In such case, 
a usage of complex automatic equalizers in the system is not economi- 
cally desirable. 

In Ref. 3, two MSE methods are discussed for the equalizer adjust- 
ment. Both methods are based on the steepest descent algorithm 
and could be easily implemented in an automatic equalizer, but not in 
a manual equalizer. 

In this paper, an MSE algorithm based on the Gauss-Seidel itera- 
tion method is described for the gain adjustment of the manual bump 
equalizers. Under the specified assumptions, this method guarantees 
the convergence of the following iterative process. If a visual display of 
the gradient of the MSE with respect to each gain setting is available, 
adjust the first gain setting until the gradient becomes zero; next, 
adjust the second gain setting until the associated gradient becomes 
zero; similarly, adjust the third gain setting and all others up to the 
last one, thus completing one iteration. As the number of iterations 
increases, the residual error in the channel will be minimized in the 
MSE sense. 

While the Gauss-Seidel iteration method may seem quite compli- 
cated, it has several distinct advantages. First of all, the Gauss- 
Seidel iteration method requires only one gradient at a time, which 
can simplify the hardware involvement, particularly for manual 
equalizers. As is shown in Section III, the number of iterations needed 
to bring the equalizers to the optimum state is not large. When the 
equalizer is composed of orthogonal networks, a single iteration is 
sufficient. Since most of the equalizer networks used in transmission 
systems are orthogonal or semi-orthogonal in nature, the number of 
iterations will usually be small. For the bump equalizers, which consist 
of semi-orthogonal terms, two to three iterations are satisfactory for 
the optimal equalization according to the error analysis given in 
Section III. This result has been verified experimentally in the field. 

II. MSE ADJUSTING ALGORITHMS FOR BUMP EQUALIZERS 

In this section, several assumptions are made before the ZF and 
MSE algorithms are presented. 
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2.1 Characterization of coaxial channel 

The channel assumed in this section is represented by an infinite 
sin x/x series on the frequency domain. (Since the transfer function of 
the Bode network is symmetric on the log / plane, where / is the 
natural frequency in hertz, the frequency used throughout is defined 
by w = log /.) 

Let M(w, t) represent the time- varying channel misalignment 
which is a real valued function of frequency in decibels. From the 
practical point of view, however, the channel can be assumed to be 
simply M (w) since the time variation is negligible during the equali- 
zation interval. Further, assume that the Fourier transform of the 
channel is limited in the time domain by a certain positive constant. 
(It should be noted that the Fourier transform of the channel does 
not result in an impulse response of the channel because the channel 
is measured in dB. However, there is an implicit dual relationship 
between the channel studied in this paper and that involving a time- 
domain equalizer, e.g., a transversal equalizer, in which a frequency 
band limitation of the channel is implied.) Hence, the channel can be 
characterized on the frequency domain by the following series : 

= £ sin[2^ - t» )] 
„=o 2rp(w - w n ) 

where C„, p, and w„ are real numbers and 

w n+ i — w n = 7T- for all n = 0, 1, • • •. 
2p 

Note that w n is equally spaced. 

Equation (1) can also be written in the following way. 

/■l =0 

M(w) = I T, C n cos[2irp(w — w„)x]dx 

J0 n=0 

— / 1 5Z C« cos(2irpio„x) cos (2irpwx) 

J0 [ n=0 

+ Y, Cn sm(2irpw n x) sin(2Trpw:r.) 



n=0 



dx 



= I {F{x) QOs(2irpwx) + H(x) sm(2irpwx)}dx, (2) 

Jo 



where 



F(x) = 21 C„cos(2Trpiu n x) and H(x) = 51 C„ sin(27rpw„a;)- 

n=0 n=0 
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Since ^ X ^ 1, eq. (2) implies that the shortest frequency domain 
ripple period found in the channel M(w) is 1/p. 

2.2 Representation ot bump equalizers 

The bump equalizer considered in this paper is a linear combination 
of adjustable-loss Bode networks. The input-output transfer function 
of an equalizer composed of N Bode networks can be represented by 






EQL(to) = L g k B k (w) (dB), (3) 



where N is the number of networks and g k and B k represent the gain 
and response respectively of the kth Bode network. 

A typical Bode network is shown in Fig. la, where the loss is con- 
trolled by the resistor R. The transfer function, B k (w), can be analyti- 
cally derived and, with a suitable flat gain amplifier, it can be expressed 
by the following equation : 

R (<n\ - [£*(! + gj) + D k J ~ D k (dB) (4) 

B k (w) - [(1 + Ek y + Dk y (dB), (4) 



where 



(w/w k ) 2 H k 

L)k = 



[(wM) 1 - i] 2 ' 



u - ^* 

"k — Y~ > 

Li k 



and 



"* - log *zm 



Since eq. (4) shows B k (w) to be a quite complicated function of w, 
one of the following assumptions is used while analyzing the equalizer 
in detail. 

Assumption 1: Let B k (w) be approximated by 

. f ir , . "1 sin[ >(w - w k )/Aw~] (AT ,, (t .s 

sine —(w — w k ) = — —, T77 \PB)' W 

I Aw K J tt(w — w k )/Aw 

Since there are N Bode networks in the equalizer, which for the 
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Fig. 1 — Adjustable Bode network. 

present analysis are spaced equally on the w-axis with interval Aw 
(see Fig. 2) , then the transfer function of equalizer can be expressed by 



N r *" 1 

^QL(w) = £ g k sine — (w - w k ) 



(dB). (6) 
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Fig. 2 — Manual equalization of L5 coaxial system. 



To permit a comparison between B k (w), as represented by eqs. (4) 
and (5), the two equations are plotted in Fig. lb. The maximum dif- 
ferences between the two best matched curves are 0.165 and 0.183 
dB when \w — w k \ S Aw and \w — w k \ > Aw, respectively. 

Assumption 2: Let B k (w) be approximated by 

. [ t, x "1 sin[?r(w — w k )/Aw] 

cosine — -(w — w k ) = — —. r-7-7- 

[ AwT J ir(w - w k )/Aw 

cos[7r(w — w k )/Aw . . 

"l -4[(u>- w k )/AwY U; 

Under the same conditions listed in Assumption 1, the transfer func- 
tion of equalizer can be expressed by 



,v 



EQL(w) = T. 9k cosine I -r-(v> ~ w k ) 



[ A^ (W " Wk) } 



(dB). 



(8) 



Expression (7) is also plotted in Fig. lb, and it can be seen that cosine 
[y{w — w k )/Aw~\ approximates quite well the actual transfer function 
of the Bode network as expressed by eq. (4). The maximum differ- 
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ence between the two best-matched curves is 0.0327 dB when \w — Wk\ 
fg Aw and 0.0404 dB for \w — Wk\ > Aw. 

2.3 Mean-squared error algorithms (Gauss-Seidel iteration method) 

The definition of optimal equalization used in this paper is the mini- 
mization of the MSE of the equalized channel as a result of adjusting 
the gain parameters, gf*. 

On a decibel scale, the equalized error will be 

E(w) = £ g k B k {w) - M(w). (9) 

Then the MSE can be represented in the frequency domain by 

MSE = f \E{w)\Hw. (10) 

J — 00 

Theorem 1 : If the equalizer described by eq. (3) is composed of linearly 
independent netivorks, then there exists a unique set of gk's which nulls 
all the gradients Gkl i.e., 

G k = d -^- =0 forallk= 1,2, ■■■, N, (11) 

ogk 

where MSE is defined in eq. (10), and the corresponding set of g^s 
results in minimum MSE. 

The proof is given in Appendix A. 

The bump equalizers considered in this paper belong to the class 
of the equalizers denned in Theorem 1. 

As derived in eq. (25) of Appendix A, the gradient vector is given by 

G = Bg - M, (12) 

where B, g, and M are the system matrix, gain vector, and correlation 
vector, respectively, and are defined as follows. Defining an inner 
product 

(A, B) = ^ A(w)B(w)div, 

J —to 

G = [G\, G2, • ■ • , Gn^\ t , 

where T indicates the transpose, 

g = [0i,02, ■ •, gN~] T , 
M - 2[<B 1 H, Mixv)), (B t (w), M(w)), ■■■, (B N {w), M(w))Y, 
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and 

\{B h B,), (B h B 2 ), ■ ■ ■ (B h B N ) 
fi = 2 (B 2 , B 1 ), (B 2 , B 2 ), • • • (B h B N ) 

(B N , Br), (B N , B 2 ), • • • (B N , B N ) 

Now the gain vector is 

g = B-i(G + M), (13) 

provided that B _1 exists. The optimum gain, g*, which results in the 
minimum MSE is obtained by solving eq. (13) with G = 0. The 
MSE algorithm given in Ref. 3 solves (13) with G = by the steepest 
descent method, which can be readily implemented in an automatic 
equalizer control circuit. For manually adjusted equalizers, however, 
the steepest descent algorithm cannot be easily implemented because 
the algorithm requires simultaneous adjustment of all the gain settings. 
A manual equalizer adjustment algorithm should have the following 
properties : 

(i) The several gain settings can be adjusted one at a time, within 

a specified sequence. 
(ii) Repeating step (i), g approaches g*. 

The converging rate of the initial g to the final g* depends on the type 
of algorithm used and the system matrix B. For the bump equalizer, 
this question is discussed in Section 3.2. 

Theorem 2 (Gauss-Seidel iteration algorithm): If the system matrix 
B in (12) has dominant diagonal elements such that 

\{Bu,B k )\ > £' \(B k ,Bi)\ (14) 

for allk=l, 2, • ••, N, 

where £' indicates the summation of all terms excluding the case j = k, 
then every g converges to the optimum gain g* = B -1 M by the following 
iteration process: 

Iteration 1: Let gf*«) indicate the kth gain at the ith iteration; thus, the 
initial gain settings are gi {0) , g 2m , ■ • •, ginoy Adjust guo) until its corre- 
sponding gradient Gi = and designate the resultant gain Qui). The gain 
settings are then gnu, g 2W , ■ ■ ■, gNm- Adjust g 2W until the gradient 
G 2 = 0, resulting in the gain settings giw, gui), 01(0), •••, 0tf(o). 
Repeating the operation for each setting results in gui), gza), Qaa), • • •, 
gN(i), and completes the first iteration. 
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Iteration 2: Adjust gnu until G\ = 0, resulting in the gain settings Qua), 
02(i), 03(1), • • ■, 0jr(i). Obtain g U2 ), g»w, • • ■, £at(2) by similar operation, 
completing the second iteration. Similarly, iterations 3, 4, • • ■ , n can be 
carried out as required. 

The proof is given in Appendix B. 

If equalizer networks satisfy the inequality (14), they are called, 
in this paper, semi-orthogonal terms. The bump equalizers denned in 
this paper satisfy the inequality (14), and hence Theorem 2 can be 
used as a manual-equalizer adjusting algorithm. To implement this 
algorithm, a visual display of each gradient is required before the cor- 
responding gain is adjusted. The following two theorems provide 
simple ways of determining the gradient. 

Theorem 3 (ZF algorithm): Let the channel be represented by eg. (1) and 
the equalizer satisfy assumption 1 . If the interval Aw between two adjacent 
Bode networks is no greater than half the shortest ripple period ( = 1 /p) 
found in the channel, i.e., 

Aw S £ , (15) 

then the optimum gain setting is obtained by repeating the Gauss-Seidel 
iteration process defined in Theorem 2 with the gradient given by 

G k = 2E(w k ), (16) 

where k — 1, 2, ■ • • , N and E(wk) is the frequency domain error value 
measured at frequency Wk, the center frequency of the kth Bode network. 

The proof is given in Appendix C. 

Thus, if signals are transmitted at a set of frequencies equal to the 
center frequencies, w k , of the Bode networks, and if the errors are 
measured at these frequencies at the receiving station, the gradients, 
Gk, can be obtained directly. Then each gain, g k , would be adjusted 
until its gradient, Gk, reduced to zero. This is the well-known "zero- 
forcing" technique used in Ref. 2; it also achieves the optimum equal- 
ized channel in the MSE sense, if the stipulated assumptions apply. 

Theorem 4 (MSE algorithm): Let the bump equalizer in this case satisfy 
assumption 2 and assume that the interval, Aw, between adjacent Bode 
networks is no greater than the shortest frequency domain ripple period in 
the channel, i.e., 

Aty < - . (17) 

V 

Then the optimum gain setting is obtained by repeating the Gauss-Seidel 
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iteration process with the gradient in this instance given by 

G k = P?( w k - ^f W ^(u;*) + i^( to* + ^ ) (18) 

fc= 1,2, ••-,#, 

where E(wk) is the frequency domain error at the center frequency of 
the kth Bode network, and E\_w k — (Aw/2)] and E[w k + (A/2)] are the 
frequency domain errors measured at lower and upper frequencies midway 
between adjacent Bode networks. Equation (18) is derived in Ref. 3. 



Proof: In this case, 
and 



(B k , B k ) = 0.75 



£' \{B kl B s )\ =0.25 



for all k, thus satisfying the inequality (14). Hence, the Gauss-Seidel 
iteration process converges to the optimum gain settings. 

To implement the MSE algorithm, a measure of the error at 2N — 1 
points in the frequency domain is required (see Fig. 2). In practice* 
the MSE technique results in better equalization than that obtained 
by the ZF method. Note that assumption 2 for the MSE algorithm 
approximates the actual equalizer more precisely than assumption 1 
does. Moreover, inequality (15) for the ZF algorithm derived in this 
section is a conservative assumption. The channel ripple period allowed 
by the MSE algorithm can be half the period assumed by the ZF 
algorithm. 

III. CONTROL OF MANUAL BUMP EQUALIZERS 

In this section, the Gauss-Seidel iteration process derived in the 
previous section is applied to the manual equalizer for optimum gain 
control. The number of iterations required to obtain acceptable gain 
settings is reflected in inequality (14). The larger the diagonal com- 
ponents [left-hand side of (14)] compared to the off-diagonal com- 
ponents [right-hand side of (14)], the fewer the iterations needed. The 
rate of convergence of the iteration is described in Section 3.2 based on 
an error-reduction analysis. It is shown that one iteration is sufficient 
to obtain the optimum gain settings for the ZF Gauss-Seidel iteration 
algorithm derived in Theorem 3. When, in the more general case, the 
channel is initially equalized by the ZF algorithm, one or two more 
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iterations are usually sufficient to achieve a practically optimum 
equalized channel for the MSE algorithm. 

3.1 Hardware realization of Gauss-Seidel iteration process 

For the L5 line equalization, an equalizer adjustment system has 
been developed for the adjustment of bump equalizers by the Gauss- 
Seidel iteration process. It is composed of a precision transmission 
measuring set, 90G oscillator — 90H detector — digital control unit 
(Ref. 4), and a hardwired, special-purpose computer which contains 
a programmed memory and an arithmetic unit called an equalizer 
adjustment unit (EAU) (see Fig. 3). Referring to Fig. 2, we assume 
that the equalizers in the receiving station are to be adjusted by the 
MSE algorithm. By selecting the particular Bode network to be 
adjusted, the EAU in the transmitting station causes the 90G oscillator 
to generate sequentially an appropriate set of three frequencies (f k i, 
fk2, and /u). The EAU in the receiving station measures the channel 
error at the same three frequencies and computes the gradient, Gk, 




Fig. 3 — Equalizer adjustment unit. 
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by the following relationship : 

G k = B k xEUki) + B k2 E(f k2 ) + B k3 E(f k3 ). 

Normally, B kl = B k3 = f and B k2 = 1. The resultant is displayed in 
a digital readout of the EAU and the operator subsequently adjusts 
g k until G k = 0. Then the next Bode network is selected and the trans- 
mitting EAU causes the generation of the required set of three fre- 
quencies, and the receiving EAU processes the received error signals 
and displays the calculated gradient. Again, the corresponding gain 
adjustment is made. When the ZF algorithm is selected, the gradient 
displayed is simply the error at the center frequency of the Bode 
network. Hence, as far as the operator is concerned, the adjustment 
procedure for the ZF and MSE algorithms is identical. 

In practice, it is found that the equalizer should be initially adjusted 
by the ZF algorithm to bring the system near the optimum state. In 
this way, any large initial gain deviations from the optimum value are 
quickly reduced to within about ±0.5 dB. Then the MSE algorithm is 
used for the "fine tuning" of the gain adjustments. Usually, one or 
two iterations with the MSE algorithm will be sufficient for the equal- 
izer to reach the optimum MSE state when starting from the ZF 
state. According to inequality (23) derived in the following section, and 
given initial gain settings within ±0.5 dB of the optimum value, the 
gain settings after two iterations are within ±0.036 dB of the ideal 
values in the worst case. The actual deviation from optimum in most 
cases will be smaller than 0.036 dB and in any case will be less than 
the inherent accuracy limitations of the 90-type transmission measuring 
sets to be used. 

3.2 Error analysis 

The Gauss-Seidel iteration provides fast convergence of initial gain 
settings to the optimum values for the bump equalizers. As derived in 
Appendix B, the Gauss-Seidel iteration can be expressed by 

g (l+1) = - IrHJgw, + L-*M, (19) 

where L is the lower triangular portion of the system matrix B in- 
cluding the diagonal and U is the upper triangular portion of B not 
including the diagonal. 

Defining an error vector after the ith iteration to be 

•<o = g* - gco, (20) 
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where g* is the optimum value, then 



•h+u =g — g(t+D 

= g* + L-»Ug (0 - L-M 



(21) 



To derive (21), the equality, L _1 M = g* + L _1 Ug*, was used. Hence, 
the magnitude of the eigenvalues of L _1 U determines the speed of 
convergence. 

If the equalizer belongs to the class defined in assumption 1, Bis 
a positive definite diagonal matrix and L -1 U is a null matrix. Hence, 
6(i + D = for all t = 0, 1, 2, • ■ •. In other words, we obtain the optimum 
gain settings by the ZF Gauss-Seidel iteration algorithm in one iter- 
ation. The same result can be obtained from eq. (19), since L^U is a 
null matrix and L _1 = B -1 . 

If the equalizer satisfies assumption 2, the MSE Gauss-Seidel 
iteration algorithm can be used. Now the system matrix is 



B = 2 



km. 



bn, 
622, 



bis, 

&23, 



biw 
bi N 



>N2, 



where 



and 



bn = 0.75 
ha = 0.125 

bn = 



1 = J 
\i ~ 3 



= 1 



if \i - j\ ^ 2 



for all if j = 1, 2, • • •, N. Splitting B into two parts, L and U, which 
are defined above, and performing some algebra, 



L- J U = 4 



0, 


6- 1 , 0, 0, 


" 1 


o, 


-6- 2 , G- 1 , 


0, 


o, 


6- 3 , -6~ 2 , 


6-S 















6 



-3 _fi-2 



6 



0, (-l)**^-*, (-1)^6-^-^, 

Hence, one can calculate the new error vector by eq. (21). After one 
iteration, the upper bound on the maximum residual error becomes 

= K6" 1 + 6" 2 + 6~ 3 + ■ 



I Ca-(I) I max 

for all j, k = 1, 2, 

|6fc(2) I mi 



) I 6j (0) I max 

= 0.26667 k< > In,.* (22) 
• •, A r . Similarly, after the second iteration, 
g 0.26667 1 e m I max = 0.07111 \e im \ m „ (23) 
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for all i, j,k= 1,2, ■■-, N. The equality in eqs. (22) and (23) will 
be obtained if and only if 

6l(0) = — ^2(0) — 63(0) = — 64(0) = ••• = |efc(0)|max- 

This in general will not be the case, and the maximum setting error 
after the second iteration usually will be less than 0.07111 times the 
maximum setting error prior to the first iteration. Consequently, if 
the gain settings are within 0.5 dB of optimum at the start (which a 
single ZF iteration will establish), the deviation from optimum settings 
in the MSE sense is within a few hundredths of a decibel after two 
additional iterations. 

IV. CONCLUSION 

This paper shows that a manual equalization process which can be 
described by a Gauss-Seidel iteration method provides optimal con- 
trol for bump equalizers in the MSE sense. Compared to the steepest 
descent method discussed in Ref. 3, the Gauss-Seidel iteration method 
can be more economically implemented for manual equalization such as 
in the L5 system. The Gauss-Seidel iteration process requires knowl- 
edge of the gradient of the MSE with respect to the gain setting for 
each Bode network to be adjusted. The ZF algorithm derived in this 
paper requires just one Gauss-Seidel iteration, but in practice the ZF 
equalized channel is not optimum and can be further improved by the 
MSE algorithm. This is because the gradient obtained by the MSE 
algorithm is more accurate than the one obtained by the ZF algorithm 
for realizable equalizer shapes. It should be noted that three tones are 
required to determine the gradient of the MSE with respect to the 
gain setting of each Bode network for the MSE algorithm, while only 
one tone is used to obtain the gradient for the ZF algorithm. The 
number of iterations that are necessary to bring the equalizer to the 
optimum state depends on how close the initial settings are to optimum. 
When the channel is initially ZF-equalized, only one or two more 
iterations are needed to optimize the channel with the MSE algo- 
rithm. This result agrees completely with experiments conducted in 
the L5 field trial. 

APPENDIX A 
Proof of Theorem 1 

Substituting eqs. (9) and (10) into eq. (11), we obtain 



G k = 2 f X B k {w) 

J —to 



3=1 



dw. (24) 
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Defining an inner product 

(A, B) = j" A(w)B{w)dw, 

0= ZG h G 2 , .»,(?»]* 

where T indicates the transpose, 



g = \jfh9h ••■,Qn~\ t 
M = 2[<fii(u>), M{w)), (B 2 (w), M(w)), ■■-, (B N (w), M(w))] T 



and 



B = 2 



(B h B 1 ), (B h B 2 ), •••, (B h B N )) 
(Bt.Br), (B 2 ,B 2 ), •••, (B 2 ,B N ) 



(B N ,Bi), (B N ,B 2 ), ■■-, (B N ,B N ) 

a simultaneous equation of the type of eq. (24) for all k from 1 to iV 
can be written as 

G = Bg - M 
or 

Bg = G + M. (25) 

Since eq. (25) is a nonhomogeneous system of N equations and G + M 
is a vector with N real-numbered components in the case considered, 
for the given G, the unknown g is uniquely obtained by 

g = B-*(G + M), (26) 

provided that B is a nonsingular matrix. 

However, if B were a singular matrix, then a linear combination of 
the columns could be made zero, i.e., 

E h k (B k , Bj) = 



for each j = 1, 2, • • •, N where hk are real nonzero numbers. 
In addition, the following relationship could also hold: 

E hMBk, Bi> + E h 2 h k (B k , Bi) + ■ ■ ■ 



+ E h N h k (B k , B N ) = 0. (27) 



But (27) can also be written as 

/ E h k B k (w) 

J-*> [ k=i 



dw = 0. 



(28) 
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Equation (28) contradicts the assumption that B k (w)'s are linearly 
independent. Hence, B is indeed a nonsingular matrix and there exists 
but one set of g k 's for which G = in eq. (26). It is yet necessary to 
prove that this stationary point is the global minimum of the MSE 
denned in eq. (10), which is established if the following relationship 
can be proved: 

r \m(w) - a Z g* k B k (w) - ^gtB k (w) Vdw 
(a J" \m(w) - Jt^tBkWVdw 

+ (3 f \m(w) - J} gTB k (w)Vdw, (29) 



where 



£ glB k (w) and £ gtB k (io) 
k=i fc-i 



indicate distinct equalizer settings, a + P = 1 and a,/3 > 0, then 
MSE is a strict convex function of gain settings g k 's and has a global 
minimum. 

Subtracting the left-hand side from the right-hand side of inequality 
(29), we obtain the following: 

-2a/? |_" I [ J] g* k B k (w)Y + [ ^ ^(w)] 2 } dw. (30) 

Since B k (w)'s are linearly independent and at least one equalizer 
setting, 

£ flEB fc (ti;) or £ ££*(«;), 

is not zero, (30) is negative. Hence, inequality (29) is correct, and the 
proof of Theorem 1 is complete. 

APPENDIX B 
Proof of Theorem 2 

The gradient of MSE with respect to the gain settings is represented 
by the following equation: 

G = Bg - M, (31) 

where B, g, and M are denned in (24). Splitting the B matrix as follows 

B = L + U, (32) 
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where L is the lower triangular portion of B including the diagonal and 
U is the remainder of B, 

G = Lg + Ug - M. (33) 

According to the iterative procedure described in the theorem, g«+u 
is obtained from g ( ,-) by setting G = 0. With the aid of eq. (33), this 
procedure can be expressed by 

= Lg (t+1) + Ug«, - M 
or 

g ( , +1) = - L-iUgc,-) + L^M. (34) 

By successive calculation, eq. (34) can be modified by 

g„+i> = C-L- J U]'+'g (0 , + t [-L-HJPL-'M, (35) 

where g(o is the initial value. 
If 

[-L-'lTl'-^CO] as ;-►«>, 

£ [-L-Wyi- 1 -* [L + U]- 1 = B- 1 . 

ft-0 

Hence, eq. (35) becomes 

g (l+ i) = + B-'M, 

which is the desired result. 

Hence, the theorem is proved if L-'U is a convergent matrix, i.e., 
eigenvalues of the matrix L-'U are all less than one in absolute value. 
However, if condition (14) is satisfied, L _1 U is a convergent matrix 
and [-L^U]' -* [0] as i -> x (see Theorems 3.3 and 3.4 in Ref. 5). 

Note: The iteration process defined by (34) is known as the Gauss- 
Seidel or, simply, the Seidel iteration. 

APPENDIX C 
Proof o/ Theorem 3 

When assumption 1 is satisfied, 

(B k , By) =1, k = j 

= 0, k^j 
for all j and k. 
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Hence, Theorem 2 is satisfied and we have to prove now 

G k = 2E{w k ). (36) 

From eq. (24), 



£ gjBj(w) - M{w) 

3=1 



dw 
w)dw 



Since 



G k = 2 f B k (w) 

J —00 

= 2 / sine — -(w-w k ) L gjBj(i 
J -oo .Am; Ji=i 

- 2 f X sine j-(w - w k ) \M(w)dw. (37) 

J sine -£-(w - to*) sine -^-(w - tty) \dw = if fc ^ j 

and 

= 1 if k = j, 

the first integration in (37) is simply 2g k . 

Substituting M(w) of (2) into (37), the second integration of (37) 
becomes 

2 f sine -£-(w - w k ) J \F(x) cos(2tpwx) 

+ H{x) sin (2irpwx) } dxdw 

= f" sine | j-(w - w k ) \ j \f(x) cos[>rp(u> - w k )x~\ 

+ h(x) sin [2irp{w — w k )x~\\dxdw, (38) 

where 

f(x) = F{x) cos (2irpw k x) + i/(x) sin (2irpw k x) 

and 

/i(rc) = H(x) cos (2irpw>*a;) - F(aj) sin (2irpw k x). 

Since 

f f h(x) sine ~-{w — iv k ) sin [2tp(w - Wfc)z]diwfa = 0, 
Eq. (38) becomes 

2 /"" sine | £-(w - w k )\ I \f(x) cos [2irp(w - w k )x]}dxdw. (39) 
Replacing w = u + w k and changing the "cos" into "exponential" 
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form, (39) becomes 
/ sine ( — u\ I { f(x)[exp(i2irpux) + exp ( — i2irpux)']} dxdu 

— I f( x ) I sine I — u) {exp(i2wpux) 

+ exp ( — i2irpux) } dudx, (40) 

where i 2 = — 1. Since ^ x ^ 1 and 2p ^ 1/Aw by assumption, 
integration of (40) is simply 

2 f 1 f(x)dx. 
Jo 

Note that the inner integration in (40) is the Fourier transformation of 
the sine function. 

Combining the results, Gk in (37) becomes 

G k = 2g k -2 f 1 f{x)dx. (41) 

Jo 

However, 

M(w k ) = \ f(x)dx and EQL(w k ) = g k . 
Jo 

Hence, (41) becomes 

G k = 2[EQL(™ fc ) - M(w t )l 

= 2E(w k ). 

This proves the theorem. 
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